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O ■ ABSTRACT 

^ \ We investigate the properties of 12 ultra-massive passively evolving early type galaxies 

^ . (ETGs) at zphot > 1-4 in the COSMOS 2 deg^ field. These 12 ETGs were selected as 

' pBzKs, have accurate 1.4 < Zphot ^ 1-7, high Sersic index profiles typical of ellipticals, 

no detection at 24^m, resulting in a complete ETG sample at M* > 2.5 x IO^^Mq 
(Chabrier IMF). Contrary to previous claims, the half light radii estimated in very 
\ high S/N imaging data from HST+ACS are found to be large for most of the sample, 

55 ■ consistent with local eUipticals. If the high redshift ETGs with M* < 2.5 x lO^M© 

are really small in size and compact as reported in previous studies, our result may 
suggest a "downsizing" scenario, whereby the most massive ETGs reach their final 
structure earlier and faster than lower mass ones. However, simulating galaxies with 
0^ ' morphological properties fixed to those of local ETGs with the same stellar mass show 

, that the few compact galaxies that we still recover in our sample can be understood in 

term of fluctuations due to noise preventing the recovery of the extended low surface 
brightness halos in the light proflle. Such halos, typical of Sersic profiles, extending 
even up to 40 kpc, are indeed seen in our sample. 
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1 INTRODUCTION 

Understanding the formation processes of galactic spheroids, 
i.e., bulges and early- type galaxies (ETG), remains a cen- 
tral theme in the context of galaxy formation and evolu- 
tion. The discovery of a widespread population of passively 
evolving ETGs at z > 1.4 proved that quenching of star 
formation in massive galaxies is well under way by z ~ 2 
jGlazebrook et alllioM : ICimatti et al.ll2004l ). 

Further studies of high redshift ETGs then revealed an 
unexpected property of a sizable fraction of them: many 
such ETGs appear to have much smaller effective radii {Rc, 
by factors ~ 2 to ~ 4) with respect to ETGs of comparable 
stellar mass (M*) in the local Universe (Daddi et al. 2005, 
DOS), a result then confirmed by several other studies (e.g. 
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McGrath et al. .2008. : Toft et al.1 12007.: ,Cimatti et al 
van Dokkum et al.l I2OO8I: iBuitragq et al.1 I2OOS 



Daraiariov et al.l |2009| : ISaracco et al.l 120091 : iRettura et al l 
200^ . ETGs with similar stellar de nsities appear to b e 
extremely rare in the local Universe l|Truiillo et al.l [20091 ) . 
although it has been argued that such compact high redshift 
ETGs hav e survived as the c ores of present-day massive 
spheroids (|Hopkins et al.|[2009l ). 

It has also been argued that the Rc — M* relation 
for ETGs would keep evolving also from 2: ~ 1 all the 
way to z ^ 0, with a steady increase of Rc for given M* 
l|van der Wei et al.| [2008: Bornardi 2009). 

The formation of very compact ETGs at high redshift 
may not be a problem: submillimetre galaxies appear to have 
comparable masses (largely in gas form) an d radii, hence 
may likely be precursors to compact ETGs jCimatti et al.l 
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I2OO8I : iTacconi et alJlioO^ . turning into them upon expul- 
sion of their residual gas. However, no generally accepted 
explanation has yet been established on how such compact 
ETGs would evolve into their present descendants, i.e., how 
they can inflate to reach 2-4 times larger effective radii. 

ETG-ETG (dry) merging soon appeare d to be a plausi- 
ble mechanism fe.g.. iKhochfar fc Silkll2006h . but most such 
events would result in an increase of both mass and radius 
such to move galaxies parallel to the local Rc — M* rela- 
tion, rather than towards it (Cimatti et al. 2008; Saracco et 
al. 2009). Moreover, major dry merging events may be too 
rare anyway, as the vast majority of high-z ETGs do not 
appear to have close ETG companions of similar brightness 
(Cimatti et al. 2008). Having excluded major mergers, ac- 
cretion of many satellites (i.e., minor dry mergers) were then 
considered, either adding an extended envelope to the com- 
pact core (Cimatti et a l. 2008), or expan ding such core by 
gravitationally heating l|Naab et al.ll2007l ). It has also been 
suggested that the expansion of high-2 ETGs would be the 
result of AG N feedback driv ing the rapid expulsion of the 
residual gas (|Fan et al.ll2008l ) , but such expansion must take 
place in at most a few dynamical times (~ 10* yr), hence 
only very few objects could be caught as already passive and 
still compact. 

On the other hand, the possibility that the apparent 
small radii may not be real was not completely excluded. 
This could be either the effect of part of the light com- 
ing from a centrally concentrated source (e.g., an AGN or 
central starburst), or to some systemati c bias not having 
been t aken in full account (DOS). Recentlv. lLa Barbera et aU 
have argued for the presence of a radial age gradient 
in local ETGs, that would result in an apparent decrease 
of Rc with increasing redshift, as the younger, more cen- 
trally concentrated population differentially brightens w.r.t. 
the older and broader stellar component. However, there ap- 
pears to be no strong colour gradients in the few cases in 
which both optica l and near-IR HS T imaging is available for 
high-2 ETGs (e.g. lToft etahlbOOTl ). or within different ACS 
bands (DOS). 

In this letter we present the results of two-dimensional 
(2D) surface brightness profile fitting of a complete sam- 
ple of 12 extremely massive ETGs (M* > 2.S x lO^A/©) 
at z> 1.4, using the HST-hACS/WFC i mages of the COS- 
MOS 2 deg^ field (|Koekemoer et al.ll2007l ). Selection of high- 
z ETGs over such a large area allows us to pick the most 
massive/luminous high-z ETGs, and tackle the size issue 
with the highest S/N ratio. 



2 DATA, SAMPLE SELECTION, AND SED 
FITTING 

A sample of extremely massive ETGs was extracted from 
the catalogue of A'-selected galaxies in the 2 deg^ COSMOS 
field (McCracken et al. 2009). From this catalogue we ex- 
tracted all the objects with Ks(Vega.) < 17.7, whose cor- 
responding BzK plot is shown in Fig. [T] From this sam- 
ple we then selected the brightes t galaxies sa t isfyin g the 
passive BzK galaxy criterion of iDaddi et al.l l|2004l ). for 
which the A's(Vega) < 17.7 limit roughly corresponds to 
M* > 2.S X 10"Mo at z ~ l.S for a Chabrier IMF. This re- 
sulted in a sample of 22 sources, having retained only galax- 
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Figure 1. The BzK two-colour diagram for all objects with 
Ks < 17.7 (Vega system) in the COSMOS field. The 12 bona- 
fide passive ETGs are shown as red circles. Blue squares with 
(z — K) < 1.6 are most likely AGN dominated galaxies, hav- 
ing point-like morphology and/or X-ray emission. The remain- 
ing blue squares with redder (2 — K) colours show sBzK galax- 
ies and pBzK contaminants rejected for being star-forming. The 
blue solid line and the red dots in the bottom show the expected 
colours for the stars (from iLeieune et al.lll997t ) , overlapping with 
the locus of real star sequences. 



ies with IRAC 3.6 detections with separations < 0".6 
from the if-band positions, a criterion required to avoid ob- 
jects affe cted by blending in the IRAC and MIPS images 
(see, e.g. iDaddi et al.ll2007l ). We notice that 4 star-forming 
BzK galaxies (sBzK) are found to K < 17.7, after exclud- 
ing likely AGNs among sBzKs, i.e., the blue squares with 
(2 — K) < 1.6 in Fig. [1] The colours of these 4 objects are 
close to the pBzK boundary and in principle some genuine 
ETG could have been scattered there just by photometric 
errors noise. However, we find that all 4 objects are well 
detected at 24 jj,m. 

Since these 22 p_B2_ftr-selected objects could include 
heavily dust-reddened star-forming galaxies at 2 < 1.4, we 
carefully checked the properties of the sample to solidly iden- 
tify the bona-fide ETGs and discard possible contaminants. 
We requested all of the following criteria for retaining ETG 
candidates: 

i) non det ection at 24 pm in the S'pi£2er-|-MIPS data 
(| Sanders et al. 2007). We used the most recent deep cata- 
logue by Aussel & Le Floc'h, reaching 5a completeness levels 
down to 80/iJy (Le Floc'h et al. 2009, submitted to ApJ). 
This corresponds to a limit of Star Formation Rate (SFR) 
< SOAf© yr~^ for 2 ~ l.S, using the models of Chary & 
Elbaz (2001); 

ii) elliptical like compact morphology, based on a visual in- 
spection of the ACS 7-band images (Fig. [2]); 

iii) multicolour SEDs best fitted with old, passively evolv- 
ing populations with no dust reddening as showing in Fig. [S] 
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SEDs were constructed using COSMOS photometry in the 
B-, i-, Z-, J, A's-band data consistently with McCracken 
et al (2009), expanded to include photometry in the first 3 
IRAC bands (from Ilbert , Salvato, & S-COSMOS collabo- 
ration, lllbert et alj (|2009l ). 

We found that 12 of the 22 pBzK galaxies do satisfy 
all these criteria and are retained as our sample of bona-fide 
2 > 1.4 ETGs. The remaining 10 galaxies that were dis- 
carded are found to fail multiple criteria in most cases: they 
tend to be 24/im detected (6/10), have irregular (or point 
like, unresolved) morphologies in the ACS imaging (7/10), 
are best fitted by star forming galaxy templates with large 
amounts of reddening (7/10), or show evidence for strong 
AGN contamination of the optical/near-IR light (3/10) as 
indicated by the 'power law' shape of the SED and X-ray 
and/or radio detection. These criteria are very efficient in se- 
lecting truly passively evolving ETGs, as confirmed by their 
spectroscopic identification in the HUDF by D05. 

Photometric redshifts and stellar mas ses were derived 
from the SEDs using the Hyperz code (|Bolzonella et al] 
l2000t ). and fitting photometric data with simple stellar pop- 
ulations (SSP) as well as with composite populations de- 
scribed by T-models (SFRoc e"'''^, with r = from 0.1 Gy r 
to oo), and using the population models of lMarastonI l|2005l ). 
This parametrisation of the star formation history may be 
inappropriate for actively star forming galaxies (especially 
at high redshifts), because it assumes that all galaxies are 
caught at their minimum SFR. However, for passively evolv- 
ing galaxies this assumption may not affect appreciably the 
redshift and mass determinations, as their SFR is indeed 
~ 0. 

The 10 galaxies that we discarded as dusty contami- 
nants are found to be at Zphot < 1.4. All 12 bona-fide ETGs 
have photometric redshifts in the range of 1.4 < Zphot < 1-7, 
are fitted by SSP models (Maraston 2005) older than ^ 1 
Gyr (i.e., with strong 4000A/Balmer breaks), and have stel- 
lar masses in the range 2.5 x 10^^ < M* < 4 x IO^^Mq, 
adopting the IChabrieJ l|2003t ) IMF, as reported in Table [T] 
Following D05 and Maraston et al. (2006) we expect that 
photometric redshifts of ETGs computed in this way should 
be accurate to within Az < 0.1. 



3 SURFACE BRIGHTNESS PROFILE 
FITTING AND SIMULATIONS 

We used GALFIT (|Peng et al.ll20o3 ) to fit 2D lSersid (|l968l ) 
profiles to the /-band (F814W filter) ACS images of the 12 
ETGs. The median integration over the COSMOS field is 
2028 s, the pixel-scale is ~ 0'.'05, and the FWHM of the PSF 
is ~ 0'.'097. GALFIT was run on s ky-subtra cted images, ob- 
tained by the SExtractor software (|Bertin fc Arnouts 1996). 
Thus, the sky was fixed to a null value in the GALFIT 
procedure, and the galaxy Sersic index (n), effective radius 
{Re), and total magnitude (mgi4) were left as free parame- 
ters. Other free parameters in the fit are the object position 
(xc, J/c), the position angle (PA) of the major axis a, and 
the axis ratio ib/a) between minor and major axis. On the 
other hand, leaving both n and the sky background as free 
parameters in GALFIT may lead to unreliable results for 
galaxies with intrinsically high values of n (see the GALFIT 



home pag43)- For each galaxy the adopted PSF was obtained 
from the closest suitable bright unsaturated star. We verified 
that swapping PSF by using different stars in the COSMOS 
field or synthetic PSFs, built by means of the Tiny TinJ3 
gives virtually identical results. When galaxy isophotes par- 
tially overlap with bright close neighbours, we fitted galaxy 
and neighbours simultaneously (this was the case for objects 
#136738, #157677, #200390, #254025). For smaller and 
fainter neighbours we used masked regions. Weight maps 
were used accounting for shot noise and sky noise. 

In Fig. [2] we show the galaxy ACS I-band images, the 
corresponding GALFIT best fit model, and the residuals for 
each object. The best- fit parameters msu, the Sersic index 
n (which measures the light central concentration), and the 
half light radius Rc (i.e. circularised radius of the ellipse 
containing half of the total luminosity of the best-fitting 
Sersic profile) are reported in Table [1] As shown in Fig. |4] 
and Table 1, only 3 of the 12 objects have very small effec- 
tive radii for their mass (0.25 < Re < 0.48 arcsec, i.e., ~ 2 
to ~ 4 kpc at 2: ~ 1.5). This qualifies them as 'compact', 
compared to the local ellipticals, whereas the remaining 9 
objects appear to be consistent with local ETGs, although 
the effective radii of some of the objects with the largest size 
(e.g., #136738, #535544, #157677) have also very large in- 
trinsic uncertainties. For all the objects but one we found 
n > 3, as expected for typical ETGs. The exception is ob- 
ject #307881, one of the smallest in the sample, which may 
consist of a bulge plus an exponential disk. On the contrary, 
the highest Sersic index (n = 7.29 ±0.58) is found for object 
#136738. A large n generally indicates an highly centrally 
concentrated light profile, compatible with such massive el- 
liptical galaxy. However, for this object the presence of some 
structure at the centre of the residual image (Fig. [2J sug- 
gests that a central unresolved point-source (probably an 
AGN) may have contributed to increase the n value. Thus, 
we fitted object #136738 by adding a central point-source, 
finding _Re=1.22 arcs (~ 10 kpc at 2; = 1.6), n ~ 4 and 
mgi4=22.28 for the Sersic component, and a magnitude of 
771814=25.88 for the point source component. We conclude 
that the possible presence of a point-source would decrease 
by almost a factor of 2 the resulting 7i value, but it would 
not dramatically affect the large Re{^ 10 kpc) found for 
this object, still consistent with the size of similarly massive 
local ETGs. 

Fig. |4] shows that our 12 ETGs span a large range in 
size, although they have very similar masses. A large scatter 
in the size distribution was also found in other studies for 
less massive galaxies. Note that the 3 smallest objects in our 
sample, with -Rc 4 kpc, are among the faintest in the full 
sample of 12 ETGs, and are all fainter than 77igi4 = 22.9. 
This leads us to wonder if the 3 small galaxies are really 
small as we have measured, or if their apparent small Re is 
due to an observational bias. Thus, what we would like to 
clarify is whether the recovered distribution of n and Rc for 
our z > 1.4 ETC sample, including the 3 apparently small 
galaxies, could still be consistent with the local distribution 
of these structural parameters. 



^ http:/ /users. ociw.edu/peng/work/galfit/TFAQ.html#sensitivity 
^ http: //www. stsci.edu /software/tinytim software pacl^age (v. 
6.3 lKrist fc Hookll2004h 
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Figure 2. ACS /-band 2.5 arcsec cutouts for the 12 bona-fide ETGs, their best-fit models and the corresponding residuals. Images are 
displayed in logarithmic scale to enhance low surface brightness details. 



In order to answer this question we have constructed 
a simulated set of galaxies all having Re = 7 kpc and n = 
6, which are the typical values for local ETGs with stellar 
masses a s large as those estimated for the 12 ETG objects at 
z > lA (jShen et a l. 2003; Caon ct al. 199;!). The GALFIT 
package was used to generate synthetic galaxies with the 
chosen light profile and having total fluxes similar to those 
of our ETG sample. The simulated galaxies were placed in 
empty sky regions of the COSMOS ACS image, close to 
our real galaxies (within ~ 1 arcmin), resulting in similar 
values of background noise. The noise and PSF images were 
constructed as for the real observed galaxies, and we tested 



that the results were not substantially altered if we used 
different PSF stars for creating and fltting the simulated 
galaxies. 

In the two panels of Fig. [5] the local average Rc — M* 
and n — Rc relations and their 1 a dispersion are compared 
with Rc, n, and masses of our 12 ETGs, along with those 
obtained for 24 simulated objects, assuming different input 
luminosities, namely mgi4,in=22, 23, 23.5, and positions on 
the frame. 

It is worth emphasising that all the 24 simulated galax- 
ies have the same intrinsic effective radius and Sersic index 
(i.e., we are neglecting here intrinsic size and n scatter). 
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Figure 3. Multicolour SED and best-fit SEDs templates (Maraston 2005, Chabrier IMF) of the 12 bona-fide ETGs. 



and yet their output 7?o,out and nout well reproduce the large 
scatter in i?c and n observed for our real sample. The spread 
in size for the simulated galaxies is a result of the spread in 
luminosity and of the different local background noise reali- 
sations, with a tendency to underestimate the size by up to a 
factor ~ 3 for the faintest objects, with m8i4,in = 23 — 23.5. 
This suggests that the apparent compactness of the 3 small- 
est objects (that are also among the faintest ones) may be 
just a consequence of their lower S/N ratio, but intrinsically 
they could not be exceptionally small in size compared to 
local ETGs of the same mass. Thus, within our sample of 12 
very high mass ETGs at < 2;phot >~ 1-6 we do not find con- 
vincing evidence for the small size syndrome affecting other 
samples of high-z ETGs. We conclude that the sample of 12 
ETGs at 2 > 1.4 shows a distribution of sizes that is fully 
consistent with those of local z = Q ETGs of the same mass. 

In order to understand these results and place them 
in the contest of previous findings we run a more general 
suite of simulations. We built 1200 synthetic galaxies with a 
much wider range of parameters: iic=0.2, 0.4, 0.8, 1.6, and 
3.2 arcsec (corresponding to ~ 2 — 30 kpc at z ~ 1.5); n=2, 4, 
6, and 8, and m8i4=20.5, 21, 22, 23, and 23.5. In particular, 
since the ACS/F814W magnitudes of the galaxies in our 
sample are in the range m8i4 ~ 22 — 23.5, this allowed us to 
explore the size-estimate dependence on the signal-to-noise 



ratio (S/N) over a range of S/N~ 50 — 200 (here as well as 
below we refer to the galaxy S/N as measured within a 1" 
diameter aperture). 

These simulations show that the size estimates is af- 
fected by a combination of three factors: i) the S/N of the 
galaxy image, ii) the intrinsic galaxy size (-Rc), and iii) the 
profile shape (n). The effects of these three factors are dis- 
cussed below: 

i) As expected, the brighter the simulated objects, i.e. 
the higher the S/N, the better the accuracy in the recovered 
parameters. As illustrated in Fig. |6l the fainter the input 
magnitude, the smaller the recovered effective radius, i.e., 
the stronger the underestimate of the size. 

ii) For simulated galaxies with parameters and S /N sim- 
ilar to those of our 12 real ETGs (msu ~ 22 - 23.5), the 
mean difference between the input and output effective radii 

turns out to be < ARc/Re > = < (-Re, in — -Ro,out)/-RG,in > — 

0.15, with a dispersion of a{ARc/ Ra) ~ 0.46, with most 
of this systematic offset and large dispersion being due to 
the largest simulated objects. For -Rc.in < 0".8 the average 
offset is very close to zero, with a small scatter. Beyond 
Ra,in = O'.'S we found that the output Ra underestimates 
the input radius and the scatter increases. 

iii) The size estimates depend on the capability of re- 
covering the Sersic index, since the best fit parameters are 
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Table 1. Characteristics of the final sample of 12 ETGs. 



#ID 


RA(J200) 


Dec(J200) 


"1814 


Jio (arcs) 


iic{kpc) 


n 




log{M*) 


age 


Av 


Ks 


93409 


10:01:27.0044 


01:39:53.1198 


22.73±0.08 


0.78±0.11 


6.70±0.92 


4.26±0.34 


1.52 


11.45 


2.40 


0.08 


17.68 


112083 


09:59:11.4551 


01:42:25.2173 


22.93±0.05 


0.44±0.04 


3.73±0.32 


3.33± 0.21 


1.61 


11.63 


2.75 


0.00 


17.51 


125656 


10:00:22.9577 


01:44:16.8353 


22.53±0.07 


0.73±0.09 


6.26±0.81 


4.37± 0.35 


1.43 


11.48 


2.75 


0.08 


17.66 


136738 


09:59:47.0280 


01:45:52.2360 


21.98±0.11 


2.02±0.45 


17.27±3.84 


7.29± 0.58 


1.59 


11.64 


2.40 


0.00 


17.35 


157677 


10:01:34.2920 


01:48:32.0022 


22.51±0.20 


3.20±1.06 


27.36±9.03 


4.46± 0.62 


1.70 


11.58 


2.50 


0.00 


17.64 


200390 


10:01:01.9974 


01:54:32.4251 


23.43±0.11 


0.68±0.14 


5.81±1.24 


4.57± 0.60 


1.61 


11.47 


2.00 


0.08 


17.68 


217431 


10:02:39.5288 


01:56:59.1162 


21.85±0.12 


2.22±0.47 


18.92±3.98 


4.81± 0.44 


1.61 


11.57 


2.40 


0.00 


17.56 


254025 


10:02:28.4912 


02:02:13.6899 


22.50±0.06 


0.68±0.07 


5.79±0.61 


4.14± 0.26 


1.71 


11.64 


2.10 


0.08 


17.54 


307881 


10:02:35.6396 


02:09:14.3640 


22.95±0.04 


0.37±0.02 


3.16±0.20 


1.92± 0.12 


1.52 


11.53 


2.50 


0.08 


17.59 


399354 


10:02:49.4055 


02:21:47.4012 


22.98±0.14 


0.75±0.20 


6.39±1.60 


4.57± 0.64 


1.70 


11.59 


2.30 


0.08 


17.63 


472257 


10:01:06.0754 


02:31:35.1966 


23.04±0.04 


0.25±0.02 


2.13±0.14 


3.63± 0.23 


1.52 


11.46 


2.40 


0.00 


17.69 


535544 


10:00:00.6592 


02:40:29.6220 


22.29±0.09 


1.06±0.16 


9.03±1.39 


4.76± 0.35 


1.41 


11.41 


2.30 


0.16 


17.54 



From left to right: #ID; identification number in the catalogue of McCracken et al (2009); RA(J2000) and Dec(J2000): 
coordinates from WIRCAM//C-band data set; mgi4; ACS F814W /-band magnitude; Re', effective radius in arcseconds and 
kpc; n: Sersic index; Zpj,: photometric redshifts, log(M*): galaxy stellar masses in Mq units; mean stellar age in Gyrs; Ay 
reddening parameter; Kg-, ii'-band WIRCAM+CHFT magnitude in Vega system). 
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Figure 4. Mass vs size for ETGs at z ^ 1.4, compared with the 
local relation (solid line) and its 1 a d ispersion (dashed lin es). 
The figure is adapted from Figure 15 of ICimatti et al.l (|2008|) by 
adding our 12 ETGs, shown as red large filled pentagons. The 
open symbols represent a compilation of th e literature res u lts fo r 
passive galaxies: big open circles are fromlCimatti et al.l (|2008| ) 
(GMASS-CDFS); op en triangles from lLonghetti et alj ( |2007| ) and 
iTruiillo et al] | |2006|) ; open square are the HUDF passive-BzK 
of Daddi et al. (2005) and Maraston et al. (2006); open dia- 
monds show the sample of van Dokkum et al. (2008) after a 
mass-rescaling of a factor ~ 1.4, to Maraston et al. (2005) mod- 
els, for homogeneity with t he other galaxies i n the figure. The 
small open circl es are from IZirm et al.l ||2007^ fFIRES-HD FS), 
Toft et al.l i200t) (F IRES- MS1054). lMcGrath et ahl l l200d) . and 
Cimatti et al.l ll2004h (K20 survey). 



corre l ated via the n — Re an d n — M* relations (|Caon et al.l 
ll993l : lD'Onofrio et al.|[l994r ). Thus, for input n > 4 the out- 
put n values tend to be underestimated up by a factor of 
~ 1.5, and mostly for galaxies with the lowest S/N. Hence, 
this underestimate of n is accompanied by a size underesti- 
mate, as illustrated in Fig (6) The recovered effective radius 




11.4 11.5 11.6 
log(MVMj 









10 






1 


, ...o- ■■■ 





1 10 
n 



Figure 5. In the top panel the stellar mass-size relation for our fi- 
nal sample of 12 ETGs (big open diamonds) is compared with the 
results obtained for 24 simulated objects (filled red symbols), all 
characterised by input Sersic parameters of nin=6, Re in ~ 7 kpc 
and magnitude mgi4 in=22 (red triangles), 23 (red diamonds), or 
23.5 (red squares). The bottom panel shows the Sersic index (n) 
- effective radius {Re) relation. The symbols are the same of the 
top panel. In both of the diagrams the l ocal M*~Re and n— He 
relations with the ir 1 a dispersions, from lshen et al] l l2003t) and 
ICaon et al. I lll993h . respectively, are represented as red solid and 
dashed lines. The black filled circle shows the input n and Re 
values common to all the 24 simulated objects. Note the large 
scatter in measured effective radius and Sersic index for the sim- 
ulated objects, well reproducing the scatter observed for the real 
galaxies. 



decreases with increasing n index (and decreasing luminos- 
ity). Similar trends are generally found for simulated galax- 
ies with i?e,in 3.5 kpc. 

These experiments confirm that for objects with large 
effective radius and Sersic index, as ETGs with masses of > 
2.5 X 10^^ M© are expected to be, with the typical S/N of our 
12 ETGs (i.e. S/N^ 50 - 100) one could stiU substantially 
underestimate n and Re- This may be due to the fact that 
the halos in high Sersic-index profiles are very extended and 
therefore have very low surface brightness, so that they can 
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Figure 6. Simulated galaxies. -Ro,out/ficin (kpc) as a function 
of input magnitude (mgi4 in=20.5, 21, 22, 23, 23.5), for a fixed 
large value of effective radius, i.e. ~ 0.8 arcsec, correspond- 

ing to -Re, in ~ 7 kpc at z ^ 1.5, and different values of input 
Sersic index, as labelled in the diagram. Worth noticing that the 
Re. out estimate decreases with decreasing luminosity and increas- 
ing Sersic index of the source. 




Figure 7. Left: the ACS I-band image of #136738 (see also 
Fig.[2j, the most massive and K-hand brightest ETG in our sam- 
ple. We have masked the central bright core and convolved the 
image with a FWHM= O'.'S Gaussian kernel in order to enhance 
low surface brightness details. A very extended low surface bright- 
ness halo is clearly detected, surrounding the galaxy. Right: the 
original ACS I-band image without smoothing. 



be easily missed in the noise at high redshift given that the 
surface brightness scales as (1 + z)"^. 

A remarkable exemplification of this situation is shown 
in Fig. [7] for the galaxy #136738, the most massive and 
K-band brightest ETG in our sample. A giant low-surface 
brightness halo is clearly detected, extending over some 4" 
in diameter around the galaxy (about 34 proper kpc sd, z = 
1.5), at a typical level of the order of 26 mag arcsec"^. Given 
the faintness of this structure, it would be obviously very 
easy to miss in the noise similar halos at a lower S/N ratio. 



4 DISCUSSION AND CONCLUSIONS 

Over the 2 square degrees COSMOS field, we have identi- 
fied and studied the 12 A'-band brightest and most massive 



(M > 2.5 X 10" Mq) ETGs at z > 1.4. Out of them, 9 
follow the same size-mass relation of local ellipticals, and 3 
appear to have smaller effective radii. However, our simu- 
lations show that their effective radius may have been un- 
derestimated by a factor of 2-3 because of their lower S/N. 
Hence, for the whole sample we find no compelling evidence 
of departures from the local size-mass relation. 

Most previous studies of high-2 ETGs were restricted to 
much smaller areas, hence they may have missed the bright- 
est/largest ETGs similar to those reported here, as their sur- 
face number density is very low, i.e., ~ 6.3±1.8 deg~^ to the 
limiting magnitude of the present sample. Using the redshift 
range 1.4 < z < 1.8 to estimate the cosmological volume, 
this corresponds to a space density of 1.4±0.4 x 10~^ Mpc"''. 
This is roughly 10% of the space density of local ETGs se- 
lected with M > 2.5 X 10" Mq. It is not clear if the different 
results of the present work compared to previous studies of 
high-2 ETGs are due to the galaxies in our sample being 
much more massive than those in most previous works. Com- 
bining the results of our 12 massive ETGs with the literature 
results -so to span a wide range of masses (2 dex)- we note 
that the departure of ETGs at 1.4 z ^ 2 from the z = 
size-mass relation seems to be differential with mass, with 
Re{z = Q)/Rb{z) being on average close to 1 for the most 
massive systems, and ~ 2 — 4 for lower masses. This is remi- 
niscent of the general downsizing scenario, whereby massive 
ETGs formed their stars and assembled their mass more 
rapidly and at higher redshifts compared to the less massive 
ones. Thus, the present results may extend the downsizing 
effect to the internal structure of ETGs, with the most mas- 
sive ETGs reaching their internal and dynamical structure 
earlier and faster than lower mass ETGs at the same red- 
shift. 

However, we also argue that bias in the galaxy pro- 
file fitting could significantly affect the size measurements, 
as low-surface brightness halos could remain undetected in 
high redshift ellipticals if observed with low S/N ratios (e.g.. 
Fig. [T]). Based on re dshifted mock images of local systems, 
iHopkins et al] l^2009^ also concluded that there is a tendency 
to underestimate the galaxy size which is more pronounced 
in the most massive galaxies, due to their large Ra and n. 
This is in agreement with the results of our simulations. 
They also stressed that elliptical galaxy profiles are not per- 
fectly "Sersic" , which would contribute to the size underes- 
timate at high redshift. 

Thanks to the exceptional brightness of our 12 ETGs, 
their S/N ratios in the single-orbit COSMOS /-band im- 
ages are among the highest reached by HST-fACS imag- 
ing ofJn£l^z_ETGSjand are comparable to those obtained 
by iDaddi etHI (|2005l ) for pBzKs in the extremely deep 
HUDF images. On average, their S/N ratio is also ~ 3 times 
higher than that of galaxies i n the GMASS stud y, which 
used GOODS-S z-band data jCimatti et al.ll2008l '). Hence, 
it is possible that also in these studies the effect of the noise 
has led to un derestimate th e size of some object by a fac- 
tor of 2 or 3. ICimatti et ahl l|2008h used simulated objects 
to check the reliability of their derived effective radii, and 
estimated an average scatter on the galaxy size measure- 
ments of ~ 20%, with no strong bias. This suggests that a 
size underestimate is unlikely to account for the small size 
obtained for the bulk of the GMASS high-z ETGs. StiU, the 
parameter space (n, Rc, S/N) explored in the simulations 
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of ICimatti et all (120081 ) may not fully include that occupied 
by the real GMASS galaxies, hence we cannot exclude that 
also in this study the size of some of the intrinsically largest 
galaxies may have been underestimated by more than a 20%. 

On the other hand, some s tudies using jJ-band 
(F160W) NICMOS-2 irn aging (e.g. iLonghetti et al] l2007l : 
Ivan Dokkum et "aP l2008l l have also reached high S /N ra- 
tios, thanks to them sampling longer rest-frame wavelengths 
where the flux of ETGs is much stronger than in the 
near-UV. Using the NICMOS Exposure Time Calculator 
(ETC), we e stimated that the S/N ratios of galax ies in 
ILonghetti eTa l. (2007) and Ivan Dokkum et al"] lj2008h were 
on average quite comparable, or sometimes even slightly 
higher than the ones obtained for our COSMOS ETC sam- 
ple. The masses of ETGs in the former study (open triangles 
in Fig. |4]) are similar to those of our galaxies, and yet the 
derived effective radii are up to 3 times smaller than those of 
local ellipticals with similar masses. Based on simulations, 
the authors argued that their typical image resolution and 
S/N led to underestimate by ~ 25% the effective radii of 
galaxies with an intrinsic n = 4 and radii in the range of 
~ 1 — 4 kpc, and took into account this result in their galaxy 
size estimate. However, they did not verify the robustness 
of their fitting procedure in recovering the light profiles of 
galaxies with i?o > 4 kpc and n > 4, typical of such ex- 
tremely massive ETGs (M > 2.5 x fO" Mq), which we 
found instead to be affected by a substantial bias. 

Similarly, Ivan Dokkum et all (|2008l ) found their ETGs 
~5 times smaller compared to local ETGs of similar mass. 
However, their stellar m asses were estimated using the stel- 
lar population models of lBruzual fc CharlotI (|2003l ) that do 
not include stars in their TP-AGB phase, and therefore stel- 
lar masses of ETGs at z > 1.4 might have bee n overesti- 
mated by a factor ~ 1.5 — 2 jMaraston et al.ll2006l ). while the 
TP-AGB effect is negligible for local ETGs. In the size-mass 
diagram, this mass re-scaling would reduce by ~ 3 times the 
average offset in size with respect to the local relation (as 
shown in Fig.|4]). Moreover, allowing for a systematic under- 
estimate of the size of the galaxies with lower S/N, similar 
to the one found by Longhetti et al. (2007) for lower-redshift 
ETGs in NIC2 images, we conclud e that some galaxy in the 
study of Ivan Dokkum et al.l l|2008l ) may actually follow the 
loca l size-mass relation. 

ISaracco et al] ('2009') collected 32 high-redshift ETGs 
with archival NICMOS data, among which 13 have given 
size comparable to the local galaxies and the others are more 
than la smaller compared to the local relation. Also in this 
case, based on the above arguments, we cannot exclude that 
the size of the galaxies with lower S/N and lower resolu- 
tion (those observed only with NICMOS-3) may have been 
underestimated. However, the authors suggested an age- 
dependent evolution, with the young ETGs (mean stellar 
age < 2 Gyr) betng more similar to their local counterparts, 
and the old ETGs being denser. Our results cannot confirm 
such claim, as our best-fit ages are all confined within a very 
narrow range (see Table[^. 

A potentially powerful test of the reliability of the size 
estimates of high-z; ETGs consists in measuring their ve- 
locity dispersion a*, which would allow us to check for the 
consistency of the three quantities entering the virial rela- 
tion: Mdyn — 5 X -RcfJ /G, with the constraint of the stellar 
mass being less than the dynamical mass. Due to the limits 



of the present-day telescopes, the low S/N of the available 
spectra of ETGs at z > 1.4 allowed only a few attempts in 
this direction, most of which based on the analysis of av- 
erage stacked spectra of sev eral objects (ICenarro fc Truiilld 
l2009l : [Cappellari et al.ll200^ ). Both these attempts are based 
on the GMASS sample of spectra l|Cimatti et al.ll2008l '), and 
within the errors found ag reement between stellar and dy- 
namical masses. However, ICappellari et al] (|2009l ) found a 
trend to underestimate the virial masses (Mdyn) on aver- 
age by ~ 30%, that may be ascribed to a size underes- 
timate by a similar quantity. Thanks to a > 7 S/N ra- 
tio of their spectra, Cappellari et al. also measured the cr^ 
from the individual spectra of two of the galaxies with the 
largest size in the GMASS sample (whose sizes are consis- 
tent with those of local ETGs). The resulting agreement 
with the masses and sizes previously estimated by Cimatti 
et al. (2008) confirmed that these two objects are truly sim- 
ilar to local ETGs, as most of the ET Gs in our sample. 
Most recently, Ivan. Dokkum et"al] (|2009l ) analysed the near- 
IR spectrum of an ETC at z = 2.2 finding a very high ve- 
locity dispersion {cr-^ ~ 510 km/s), consistent with the small 
size measured for this object, which then suggests a strong 
structural/dynamical evolution from z ~ 2 to for similarly 
massive ETGs. Yet, the relatively low S/N of the near-IR 
spectrum implies a fairly large uncertainty in the measured 
(T^. Recently some among us were able to obtain a good S/N 
near-IR spectrum of one of our 12 ETGs (object # 254025) 
with the MOIRCS spectrograph at the SUBARU telescope, 
and a measurement of cr* will be attempted (M. Onodera et 
al. in preparation). 

We conclude that the fraction of high-redshift ETGs 
which are genuinely undersized with respect to local ETGs 
is still poorly known. Several massive ETGs at high redshift 
(certainly the majority in our sample) are as large as their 
local counterparts, and therefore recipes aimed at inflating 
putative small ETGs should avoid to create oversize ETGs 
at low redshifts out of those having already reached their 
final size at 2: > 1.4. We also hint at a possible downsizing 
effect not only on stellar ages and mass assembly, but also 
on the structural properties of these galaxies, with the most 
massive ETGs being in all respects the first to reach their 
final configuration. 
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